Relaxation, pre-thermalization and diffusion in a noisy Quantum Ising Chain. 
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We study the dynamics of thermalization resulting from a time-dependent noise in a Quantum 
Ising Chain subject to a sudden quench of the transverse magnetic field. For weak noise the dynamics 
shows a pre-thermalized state at intermediate time scales, eventually drifting towards an asymptotic 
infinite temperature steady state characterized by diffusive behavior. By computing analytically the 
density of kinks, as well as the transverse and longitudinal magnetic field correlators, we characterize 
these two regimes, their observability and their signatures in the various physical quantities. 
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The dynamics of relaxation towards thermal equilib- 
rium has been one of the recurrent themes of theoreti- 
cal physics in the past decades [1] |2]. The problem is 
of crucial importance in many contexts, ranging from 
condensed matter physics to cosmology: if we think of 
injecting suddenly, e.g. by an abrupt change of one of 
its parameters (a quantum quench)^ a finite amount of 
energy in an otherwise closed many-body system, un- 
der which conditions will the system reach a thermal 
steady state ? And how is the steady state going to 
be attained ? The first question has been thoroughly 
addressed in the literature [1-4 : on one hand it is nat- 
ural to expect that scattering processes will in the long 
run lead to an ergodic, thermal redistribution of energy 
among the elementary degrees of freedom [3^ • An excep- 
tion are however integrable systems, where multi-particle 
scattering processes are highly constrained as a result of 
conservation laws [4 . As recently observed in experi- 
ments with quasi-ld Bose gases, thermalization in the 
usual sense will not occur [5 and the asymptotic state 
eventually attained by the system in the thermodynamic 
limit is expected to be described by an effective General- 
ized Gibbs ensemble (GGE) accounting for all conserved 
quantities [4[[6]. 

While in the past few years a great deal of attention 
has been paid to the description of the asymptotic steady 
state, much less is known about the dynamics of equili- 
bration of both thermally isolated and open quantum 
many-body systems. A common feature in both cases 
is the expectation of a the dynamics towards the steady 
state is characterized by various stages. In the case of 
open quantum systems, such as spin chains coupled to 
classical and quantum uniform noise [71 18] or to a bosonic 
bath [Hi [To] , these are expected to be driven by the in- 
terplay between many-body interactions, dephasing and 
dissipation associated to the external environment [TT] . 
For example, the thermalization of the system follow- 
ing a rapid quench of the bath temperature is driven 
by the spread of thermal correlations at a velocity set 
by the temperature ^Oj. For closed systems recent 
studies of the dynamics of quantum field theories [T2] 
suggest that first the system decays to a so-called pre- 
thermalized state, where the expectation value of cer- 



tain macroscopic observables is to a good approximation 
thermal, while the distribution function of the elemen- 
tary degrees of freedom is not [12] [13]. Only at a second 
stage, when energy is efficiently redistributed by scatter- 
ing processes, real thermalization occurs. Signatures of 
these crossovers have been investigated theoretically in a 
variety of systems (e.g. Fermi-Hubbard models [13f[T5] 
or Spinor Condensates [16 ), and prethermalization has 
been observed experimentally in split one dimensional 
condensates [171 [IH]- While it is evident that the dy- 
namics of thermalization could in general display vari- 
ous crossovers, including a pre-thermalized plateau, it is 
not clear from the outset whether this is a general phe- 
nomenon for both open and closed systems, what are the 
conditions for its observability on the system at hand 
and what are going to be the signatures in the various 
observables. 

In this Letter we address these issues by considering 
the dynamics of thermalization of a prototypical weakly 
perturbed integrable system, a Quantum Ising Chain 
subject to time-dependent noise. The perturbation con- 
sidered in this work does not conserve the energy, hence 
the dynamics of quantities averaged over the noise re- 
sembles that of an open quantum system. Studying an- 
alytically the dynamics of all essential observables and 
correlation functions following a sudden quench of the 
transverse magnetic field in the Quantum Ising Chain, 
we will discuss the nature of the thermalization dynam- 
ics and show that it is characterized by a crossover be- 
tween pre-thermalized and thermalized regimes. In par- 
ticular, we will show that pre-thermalization originates 
from the spreading of quantum and thermal correlations 
at different velocities. This effect is clearly observable for 
the transverse magnetization, where it leads to a neat 
crossover towards a long time diffusive behavior of the 
correlators. On the other hand, it leaves much weaker 
signatures on the correlation functions of the order pa- 
rameter, which has always a "thermal" form ^^. Finally 
we thoroughly discuss which features of the dynamics of 
this open system are expected to be relevant for the more 
complex problem of thermalization resulting from inte- 
grability breaking in a thermally isolated system. 

Before entering technical details, let us summarize the 



main qualitative picture emerging from our analysis. We 
are going to study a weakly perturbed Quantum Ising 
chain, characterized by the hamiltononian H = Hq + F, 
where Hq 



Ho 
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describes the Integrable Quantum Ising chain [20 . Here 
cr^'^ are the longitudinal and transverse spin operators at 
site i and g is the strenght of the transverse field, while 
V = J^^ Sg{t)af is a weak time-dependent white noise, 
with zero average and a strength characterized by the 
parameter F, 



{5g{t)5g{t')) = -5{t-t'). 



(2) 



The Quantum Ising chain is among the simplest, yet non- 
trivial integrable many-body system, whose static [20] 
and dynamic properties [19l |2T] are to great extent 
known. It is characterized by two dual gapped phases, 
quantum paramagnetic (^ > 1) and ferromagnetic {g < 
1) separated by a quantum critical point (^ = 1) where 
the gap A = \g — l\ closes. For a quench of the transverse 
field all essential correlation functions have been studied 
extensively ^9\ |2T]. In the following, we will consider 
the dynamics of the noisy Ising chain following a quench 
protocol: for times t < the system is assumed to be in 
the ground state of Hq with g = go and Sg{t) = 0, and at 
time t = the noise is turned on together with a global 
quench of the transverse field go ^ g. We will present re- 
sults for quenches within the paramagnetic phase - other 
types of quenches will be discussed elsewhere [22^ . 

A qualitative picture of the mechanism of pre- 
thermalization can be obtained by looking at the sepa- 
ration of time scales associated to three distinct phys- 
ical effects [23]. First of all, the coherent superpo- 
sition of modes with different frequencies leads to a 
first, power law decay of physical quantities towards the 
pre-thermalized (yet non-thermal) state described by a 
GGE [12] [151 [24]. This effect is similar to inhomogeneous 
broadening. This first decay is abruptly accelerated on 
time scales of order of 1/F by the intervention of noise in- 
duced dephasing. Finally, on much longer times scales the 
occupation of quasi-particles evolves towards its thermal 
value (in our case corresponding to infinite temperature). 
All of this is clearly observed in the time evolution of the 
density of kinks nkink(^) = (X]i(l~<^f <^f+i)/2), which can 
be written as Ukinkif) = Udriftif) + An(t); here, Udrift 
describes the heating of the system towards the state 
of infinite temperature, and for finite F and long times 
is given by nkink — \ — , J^^/r^ i while An(t) describes 
the sharp oscillations related to the quantum evolution 
induced by the quench (see Fig. [l] and discussion after 
Eq.(16)). 

This sequence of crossovers can be further analyzed 
by studying the spreading of correlations pa{r^t) = 
{o-f{t)af_^^{t)) in the transverse {a = z) and longitu- 
dinal {a = x) directions. A pre-thermalized plateau 




FIG. 1: [Color online ] The density of kinks vs. time for 
a quench with F = 0.01, A = 3, Aq = 0.1 (blue full line). 
While the red dotted curve is the value attained by rikink 
without perturbation (the red horizontal straight line sets the 
asympthotically value predicted by GGE), the full blue line 
shows the time evolution with noise that first appears to sat- 
urate at the GGE value, but later runs away towards the 
infinite temperature state (the blue oblique straight line is 
ndrift)- [Inset ] Ankink vs. time. The dephasing of the quan- 
tum oscillations is accelerated in the presence of the noise 
(blue full line) in comparison with the case without pertur- 
bation (red dotted line). 



emerges when quantum and thermal correlations prop- 
agate at sufficiently different velocities: this is manifest 
in the various regimes of the transverse correlation func- 
tions. The first crossover is observed for Ft <C 1 



p""{r,t)-rt<^i < 



'2;^exp[-2Aor] r>t 
,2^exp[-r/^^] r<Ct 



(3) 



where ^z is the correlation lenght associated to the quan- 
tum quench protocol For larger times, namely Ft :^ 1, 
the noise becomes relevant and the second crossover, be- 
tween exponential and diffusive behavior of the correla- 
tor, emerges 
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with 7 = F/A. The correlation length, which is dictated 
at t <C r by the initial gap Aq, crosses over as r :^ t to an 
intermediate asymptotic form depending by the masses 
of the quench, through the function ^^ = (,z{go^9i) [21] . 
Though admittedly thermal behavior is stronger in the 
longitudinal correlators (see below), the state reached in 
the intermediate plateau is the one the system would 
have reached in the absence of noise, hence in a sense the 
pre-thermalized state of the weakly perturbed integrable 
model [15] |18] . The first crossover is driven by the light- 
cone effect [25l [26] : this intermediate plateau persists 



until a second front, propagating at speed lower by a fac- 
tor 7 and carrying thermal correlations, passes through 
r. At this point the correlator crosses over to a diffusive 
time-dependent form, consistent with thermal Glauber 
dynamics [27 , indicating the continuous heating of the 
system towards the infinite temperature state. A similar 
behaviour emerges also for quenches across the critical 
point or starting at the critical point. In the absence of 
the quench, the correlator doesn't show the intermediate 
prethermalized regime and the equilibrium correlator - 
g-2Aor_ (3j.Qgggg QYQY directly to the diffusive behaviour 



equilibrium state with gap Aq one obtains 
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where 1/m = (1/C + 2(1^5^) ^^^ ^/^ = ^°s(^ + ^o^" 
Since the same exponential form persists the spreading 
of quantum and thermal correlations will not result in a 
diffusive form, but rather modify just the specifics of the 
correlation length which at later times shrinks as l/Tt. 
The different signatures observed in the transverse and 
longitudinal magnetization are consistent with analogous 
effects observed elsewhere [19l [30] . 

Let us now present some of the technical details of 
our analysis. Our main task will be to write a closed 
kinetic equation for the quasi-particles of the Quantum 
Ising chain subject to the infiuence of noise. Quasi- 
particles can be introduced in terms of Jordan- Wigner 
fermions |20|, q, through the relation af = 1 — 2cJq and 
(jf = —Ylj^^ii^ — 2cjCj){ci + cj). The Hamiltonian takes 
in Fourier space, Ck = XI .• CjC^^^ , the simple form 



H 



k>0 



where 



Hk = Hi + 8g{t)cj, 



(6) 



(7) 



FIG. 2: The spreading of quantum and thermal correlations 
in the noisy Quantum Ising Model: the transverse field corre- 
lator has a first crossover when ballistic quasi-particles, carry- 
ing quantum correlations, propagate at the distance r. Ther- 
mal correlations propagate at a second stage, leading to a 
crossover to a diffusive form, consistent with thermal Glauber 
dynamics. 



Thermal behavior in the intermediate plateau is observed 
more clearly in the correlation function of the longitudi- 
nal magnetization. For a quench without dissipation this 
correlator has in general an exponential form p^^ (r, t) ~ 
exp[— r/^], with a correlation length ^ dictated by the 
non-thermal distribution function of quasi-particles and 
predicted by the Generalized Gibbs ensemble [28 . For 
small quenches, however, it turns out that ^ can be ef- 
ficiently parameterized in terms of an effective tempera- 
ture Teff [H] in the form ^ ^ ^J^l /{2T^^/\) exp[A/Teff]. 
The emergence of this thermal behavior (not to be con- 
fused with thermalization) in a system with a non- 
thermal distribution of quasi-particles is strikingly simi- 
lar to the phenomenology of pre-thermalization observed 
experimentally in split condensates [17 (the two effects 
turn out to have a common physical origin [29] ) . Turning 
on the noise, the signatures of the crossover observed for 
the transverse magnetization are expected in this case to 
be different. Indeed, computing the correlation function 
following a switching on of the noise starting with the 



{g — cos k)az — {smk)(jy^ ifjk is the Nambu 
and (jn,, (jz are the Pauli matrices in the 



Hi 

spinor (/^y, cn^ ^y, 
2x2 Nambu space. Without noise, the diagonal form 
H = T.k>o^killlk - 1-kllk)^ with energies Ek = 
y {g — COS /c)^ + sin^ /c, is achieved after a Bogoliubov ro- 
tation Ck = Uk{g)jk - ivk{gW_k ^^^ ^-/c = Uk{g)jlk - 
ivk{g)'^k] the coefficients are given by 

Uk{g) = cos{Ok{g)) Vk{g) = sm{Ok{g)) (8) 

where tan(26>/c(^)) = sin(A:)/(^ — cos(A:)). Notice that the 
gap in the spectrum is A = |^ — 1| [20] . 

In order to describe the dynamics of the Ising model 
under the effect of the noise we will now derive a set of 
kinetic equations for the lesser Green function [31] 

G<{t,t') = \G<{t,t')\ . . = i{4^j{t')iJk,i{t)), (9) 

in terms of which we will express all physical observables 
of the model. Our kinetic equation is encoded in the 
Dyson equation for the contour ordered Green's function 



GC /OC j^ /OC V^C /^C 



(10) 



where Gg / i^ the unperturbed Green function and H^ ^/ 
is the self energy; in right hand side it is understood 
a convolution product, all the quantities are evaluated 
along the Keldysh contour. In the following we will take 



all self-energies within the self- consistent Born approx- 
imation [9], controlled by the small parameter 7. Spe- 
cializing now the Dyson equation to the lesser and ad- 
vanced Green's function and defining the equal time ma- 



trix pt = 
ulations 



-iGf^^ we obtain with simple algebraic manip- 



dtpk = -i[Hl,pk] + -{(jpkcr- pk), 



(11) 



coherences in turn decay exponentially fast: focusing on 
/c !:^ 0, and using 7 <C 1, we indeed find 
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the second term in the right hand side takes into account 
the effect of the noise, and a = cos20kCFz + ^vii20}.cFy 
[32] . Here pk are expressed in the basis of the fermions 
diagonalizing Hq 



Pk 



\{i-kik) ii-kilk)^ 



(12) 



where {Tklk) are the populations of levels of momentum k 
and (7^7_;.) the coherences. While this kinetic equation 
is analogous to a Lindblad master equation for a two 
level system, the Keldysh diagrammatic technique has 
the advantage of clarifying the approximations made in 
a language appropriate, and easily generalized to other 
many-body problems pT] . Parameterizing the matrix pk 
in the form p/. = ^l + Sfk^z + Xk^x + VkO'y, we obtain a 
set of three coupled equations. 



dtXk 
dtVk 



-rsm^20kSfk 

-Txk - 2Ekyk, 



-yk sin 40k, 



(13) 
(14) 



= - sinAOkSfk + 2EkXk - Tcos^ 2(9/e^fe.(15) 



The initial time condition depends of course on the phys- 
ical situation to be studied: for a quench from ^0 > 1 
to ^ > 1, one obtains pk = 1/2 + (cos(2Aa/c)/2) a^ + 
(sin(2Aa/c)/2) (Jy (where Aak = Ok{g) - Ok{go)). These 
equations can be readily solved exactly and the resulting 
long time dynamics turns out to be dominated by the 
modes close to /c = 0, ±7r This can be clearly evinced by 
neglecting the terms coupling 5fk — Vk in Eq.([l3|-([l5|, 
in the limit 7 <C 1. In that case Eq. (13) gives Sf{t) = 
(sin^(Aa/e) - l/2)e-^^^"' ^^^S where Aak = Oi - Oq. In 
this expression we can now see that the relaxation rates 
tend to vanish close to the band edges (A: = ± tt), while 
most of the modes relax fast to their thermal occupa- 
tion (n/c c^ 1/2) on time scales of the order of l/F. The 



where we have defined p_ = (Aq — A)/Ao [3T, a c^T -\- 
0{k'^) and P :^ 2A + 0(/c^). A similar analysis expansion 
is of course possible close to k = ±7r. 

Using the above expressions for the populations and 
the coherences, the computation of the number of kinks 
and of the correlators follow the standard tecniques used 
for the Quantum Ising Model [H |20l [25l [211 [MM] 
and the results shown in the first part of the paper 
comes after ordinary algebraic manipulations, that we 
will report elsewhere [22^. For instance, in the number 
of kinks, Ukink = ridrift{t) + An(t), where Udriftit) = 



E 



1 

k>0 2 



SfkCos2Aal and An{t) = ^j^^Qsin2Aalyk, 



the first term describes the drift towards the infinite tem- 
perature state and the second one, due to the coherences, 
describes oscillations suppressed by the noise on a time 
scale ~ 1/r (see inset of Fig. 1). 

In conclusion, we have discussed the dynamics of ther- 
malization in a weakly perturbed integrable model, a 
noisy Quantum Ising chain. Though the results presented 
in this work pertain to a noisy, open quantum system, 
we believe that some of the features discussed, such as 
crossovers between prethermalized and thermal regimes 
driven by the spreading of fronts at different velocities, 
will be observed also in the dynamics of thermalization 
of weakly non-integrable models. Similarly, we do expect 
different signatures of thermalization to be observed con- 
sistently in the transverse and longitudinal correlators 
[T9[[30l. At the present, however, it is not clear whether 
a weak breaking of integrability in the Ising model will 
generically lead to diffusive behaviour of the transverse 
magnetization correlator: this issue remains to be ad- 
dressed by future studies. 
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